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Abstract. The nearby red supergiant (RSG) Betelgeuse has a complex circumstehar medium out 
to at least 0.5 parsecs from its surface, shaped by its mass-loss history within the past ~ 0.1 Myr, 
its environment, and its motion through the interstellar medium (ISM). In principle its mass-loss 
history can be constrained by comparing hydrodynamic models with observations. Observations 
and numerical simulations indicate that Betelgeuse has a very young bow shock, hence the star 
may have only recently become a RSG. To test this possibility we calculated a stellar evolution 
model for a single star with properties consistent with Betelgeuse. We incorporated the resulting 
evolving stellar wind into 2D hydrodynamic simulations to model a runaway blue supergiant 
(BSG) undergoing the transition to a RSG near the end of its life. The collapsing BSG wind 
bubble induces a bow shock-shaped inner shell which at least superficially resembles Betelgeuse's 
bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock 
generated by the now defunct BSG wind's interaction with the ISM. We investigate whether 
this outer shell could explain the bar feature located (at least in projection) just in front of 
Betelgeuse's bow shock. 

1 Introduction 

Betelgeuse (a Orionis) is one of the brightest stars on the sky and one of the nearest red supergiants 
(RSG), but its circumstellar medium (CSM) has only begun to reveal its secrets in the past 20 years. 
Part of the reason for this is that optical and near-infrared observations are very difficult because the star 
is very bright at these wavelengths, and so most detections of the CSM are from space-based far-infrared 
and ultra-violet observatories. IRAS observations (Noriega-Crespo et al. I1997|) of Betelgeuse detected a 
bow shock and a linear "bar" feature upstream from the bow shock and perpendicular to the star's proper 
motion. Later AKARI observations (Ueta et al. I2008P confirmed these features with higher resolution 
observations of part of the bow shock. Herschel observations (Cox et al. 120121 Decin et al. I2012p have 
dramatically better spatial resolution than previous data and resolve the bow shock into a number of 
layers. Herschel's better wavelength coverage also allows more detailed dust modelling, providing mass 
estimates for Betelgeuse's circumstellar structures (Decin et al. I2012p . A thin shell of far-ultraviolet 
emission was found by Le Bertre et al. (|2012p at the same position as the IR bow shock in archival 
GALEX data. They also presented VLA 21cm observations with evidence of some atomic hydrogen near 
the bow shock, and stronger emission from an inner shell three times closer to Betelgeuse. 

Mohamed, Mackey, & Langer (|2012p presented the first 3D hydrodynamical simulations of Betelgeuse's 
bow shock, assuming a constant wind. They found that for a wide range of parameters the bow shock 
should have a mass M k, 0.05 — 0.15 M©, more than lOx larger than the mass M w 0.0033 Mq they 
estimated from AKARI data. This result was confirmed in the 2D bow shock simulations of Decin et 
al. (j2012p . and the discrepancy is even more severe when one considers the best mass estimate from 
Herschel of M w 0.0024 Mg (Decin et al. [20T2)) . 

Mohamed et al. (|2012p concluded that the bow shock must be very young (< 30 kyr) and not yet in 
a steady state. This means that either the stellar wind or the ISM properties have changed dramatically 
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Fig. 1. The wind mass- loss rate, M in 10 ^ Mq yr ^, wind velocity, in kms ^, and M/v^ (proportional to 
wind density, arbitrary units) are shown. Open circles identify the times of the four snapshots in Fig. (2] 



on this timescale. They further suggested that if Betelgeusc had recently evolved from a blue supergiant 
(BSG) to a RSG, then the upstream bar could be a remnant bow shock from an earlier main sequence 
or BSG wind, currently being overtaken by the star. In Mackey et al. (|2012p we tested this idea by 
simulating the GSM produced by a runaway 15 Mq star as it evolves from a BSG to a RSG at the end of 
core helium-burning. The results of Mackey et al. ()2012p are summarised in this paper, and are further 
discussed in the context of the new results from Decin et al. (|2012p . Section [5] describes the stellar 
evolution model and simulation setup, section |3] summarises our results, and section |4] presents some 
discussion and conclusions. 



2 Description of Model 

A 15 Mq solar metallicity star is modelled travelling at = 50 kms^^ through the diffuse ISM with hy- 
drogen number density nn = 0.2 cm~^ and temperature T = 7700 K. Stellar evolution is calculated using 
the Yoon & Langer (|2005p stellar evolution code including mass loss but no convective core overshooting. 
The tip of the RSG branch in this model has properties consistent with Betelgeuse (Neilson et al. 120111 
and references therein). Wind velocities are approximated using the Eldridge et al. (|2006p prescription, 
= Pw{T)v1^^. To approximate the slow wind of Betelgeuse we use jS^, = 0.04 for T < 3600 K (instead 
of the recommended value of 0.125). After the main sequence the star becomes a RSG; it then embarks 
on an extended blue loop for most of the He-burning phase; and finally evolves again to a RSG for the 
last K, 40 kyr of its life. The wind velocity, mass- loss rate M, and wind density (AI/vw), are plotted in 
Fig. [T] for the last w 75 kyr of evolution, corresponding to the end of the blue loop and evolution along 
the RSG branch. 

The circumstellar medium is calculated by solving the inviscid Euler equations for an ideal gas on a 
uniform axisymmetric grid in (z, R) in the rest frame of the star. A second-order-accurate integration 
scheme is used, and radiative cooling is included as a source term to the energy equation using the 
coUisional ionisation equilibrium cooling curve of Wiersma et al. (|2009p . The code is described in more 
detail in Mackey et al. (|2012p and Mackey & Lim (|2010p . The stellar wind boundary condition is 
implemented by imposing a freely expanding wind within a 20 grid-zone radius of the origin, updating 
the wind parameters from the stellar evolution model every timestep. 
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Fig. 2. Log of hydrogen number density nu (upper half-plane, in cm^"^) and log of temperature (lower half-plane, 
in Kelvin) plotted on logarithmic scales at times t = 11.85, 11.875, 11.8875, and 11.9073 Myr from left to right, 
respectively. Velocity streamlines are also plotted in the upper half-plane. Only part of the simulation domain is 
shown. The snapshot times correspond to the circled times in Fig. [T] 



3 Results 

Snapshots from the simulation are shown in Fig. [2] at the four evolutionary times marked by circles in 
Fig. [TJ Log of riH is shown in the upper half-plane together with velocity streamlines traced from the 
right-hand boundary. The lower half-plane shows log of T, where the rapidly decelerating wind is seen 
by the decreasing temperature in the shocked-wind region from T > 10^ K initially to T < 10"* K in the 
last snapshot. 

The first snaphot shows the CSM at the end of the BSG evolutionary phase, with a single bow shock 
comprising (from small radii to large): the freely expanding wind out to r w 0.5 pc; the reverse shock; 
the hot shocked wind bubble at lower density than the ambient ISM; the contact discontinuity distorted 
by Kelvin-Helmholz instability; the cooled and shocked ISM shell; and finally the forward shock at about 
r = 0.95 pc along the symmetry axis. 

In the second snapshot the star has become a RSG and hence the hot bubble has lost its energy 
source, so it is cooling and expanding back towards the star. The streamlines show that the contact 
discontinuity no longer has strong shear across it because it has also begun receding to the star. 

By t ~ 11.8875 Myr in the third snapshot, the RSG wind has begun interacting with the collapsing 
hot bubble, generating an inner shell around the upstream hemisphere of the RSG wind that is shaped 
like a bow shock. The final snapshot shows the CSM at the prc-supcrnova stage. Here the hot gas 
from the BSG wind bubble is all downstream from the star, the RSG wind is shock-confined (but still 
advancing) on all sides, and the upstream shock is now advancing into undisturbed ISM. 

In Mackey et al. (|2012p we suggested that the inner shocked-shell in the third panel of Fig. [5] could 
correspond to the bow shock seen at r « 0.35 pc around Betelgeuse, and the receding BSG bow shock is 
a plausible candidate structure for the upstream bar at r w 0.45 pc. The masses of the unshocked RSG 
wind, of the shocked ISM from the relic BSG bow shock, and of the inner shocked-shell are plotted as 
a function of time in Fig. [3] during the time when the inner shocked-shell exists. The RSG wind mass 
increases almost linearly to M « 0.05 Mq over this timescale, and at the same time the mass of shocked 
ISM from the BSG bow shock decreases to about the same value as it is swept downstream (only shocked 
ISM upstream from the star is counted). The inner shocked-shell is much smaller than both of these, 
growing from zero to M « 0.002 Mq at 11.891 Myr, when it reaches the upstream contact discontinuity 
and begins to merge with it. 
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Fig. 3. Masses of CSM structures as the simulation evolves. The green solid curve shows the total mass in the 
BSG bow shock upstream from the star. The blue solid line shows the total mass of freely-expanding wind within 
0.3 pc of the star. The black solid and dotted line shows the mass of the inner shocked-shell as a function of time 
in the region upstream from the star. At t ~ 11.891 Myr the inner shocked-shell begins to merge with the contact 
discontinuity from the outer shell; this merged overdensity is shown with the dotted line. 



4 Discussion 



4.1 Comparison with recent observations 

Decin et al. (|2012p have carefully analysed the Herschel Betelgeuse observations with detailed dust 
modelling, to estimate the masses contained in the bow shock and bar. We compare here with their best 
estimates for the masses, but note that the uncertainties in the estimates are at the ordcr-of-magnitude 
level. They obtain a bow shock total mass of M w 0.0024 M©, in good agreement with the estimate of 
M w 0.0033 Mq by Mohamed et al. (12012^ from AKARI data. This is marginally larger than the mass 
of the inner shocked-shell in our simulation, but given the uncertainties this is very good agreement. 

If the linear bar is composed of ISM gas (which it is in our model), with ISM dust content, Decin et 
al. p012p calculate its mass to be M k, 0.01 Mq. This is 5 — 7x smaller than the remnant BSG bow 
shock in our simulations (see Fig.[3|), but in our model only the part of this mass closest to the star would 
be emitting brightly in the far-infrared. The main argument against our model from the Herschel data is 
that the bar seems to have no curvature, whereas our simulation shows some curvature. The rotational 
symmetry of the simulations presented here significantly limits the bow shock shape however; in particular 
it can only be perfectly linear if it is aligned exactly perpendicular to the motion. Three-dimensional 
simulations that allow the bow shock to move in all directions are needed to determine whether such a 
linear feature would naturally arise from a circumstellar bow shock. 

The mass of the RSG wind was measured with 21cm neutral H observations by Le Bertre et al. (|2012p . 
and they obtained M w 0.068 M©. With a different neutral H dataset, Decin et al. (|2012p derive a similar 
mass of M « 0.02 — 0.07 M©. These measurements are within a factor of 2 of the RSG wind mass in 
our simulations, although it is something of a mystery that Le Bertre et al. p012p see evidence that the 
neutral H is confined to a radius r < 0.12 pc. Nothing in our simulations provides such a confinement. 

The GALEX far-ultraviolet emission (Le Bertre et aZ. 12012^ at the same position as the IR bow shock 
may further support our model. If this emission is excited by hot electrons, the remnant BSG hot bubble 
is a natural place to find them, and they would interact directly with the outer edge of the expanding 
RSG wind. Further work is required in order to make a quantitative prediction of the far-ultraviolet 
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emission from our simulations. 
4.2 Summary 

We have presented a new model for the CSM around Betelgeuse, in which the star has recently evolved 
from a BSG to a RSG and is now within 20 kyr of the end of its life. The evolving stellar wind together with 
the stellar motion produce a complex CSM with multiple shells and even regions of hot gas with T > 10^ 
K within 1 pc of the star (especially downstream) . The inner shocked-shell around the expanding RSG 
wind has a comparable mass and radius to the observed bow shock around Betelgeuse, and is 20 — 50 x 
less massive than an equilibrium RSG bow shock of a similar size. 

The remnant of the BSG bow shock forms a large structure upstream from the inner shocked-shell, 
with a much larger radius of curvature because it originally had a much larger standoff distance from the 
star. Its location makes it a plausible candidate for the bar upstream from Betelgeuse's bow shock, and 
its mass seems to have the right order-of- magnitude (although this is very uncertain). Its curvature does 
not agree with the observed bar, which has no discernable curvature, but this is not robustly predicted 
by our two-dimensional simulations. 

The model we present therefore has some encouraging agreement with observations, in particular in 
explaining the very low mass of the bow shock, which is difficult to explain in any other way. On the other 
hand some questions remain, notably the apparent confinement of the RSG wind to r < 0.12 pc from 
neutral H observations, and also the detailed structure of the linear bar. Further observational study of 
the CSM is therefore required, in particular data that measure the gas as well as the dust, because the 
dust-to-gas mass conversion is very uncertain and prevents us placing strong constraints on our model. 
Spectral data that constrain the gas velocity and its chemical state would also be very useful, although 
the star is so bright that this is difficult at optical wavelengths. 
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